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Abstract 

Integral transformations of the QCD invariant (running) coupling and of some 
related objects are discussed. Special attention is paid to the Fourier transformation, 
that is to transition from the space-time to the energy-momentum representation. 

The conclusion is that the condition of possibility of such a transition provides us 
with one more argument against the real existence of unphysical singularities observed 
in the perturbative QCD. 

The second conclusion relates to the way of "translation" of some singular long- 
range asymptotic behaviors to the infrared momentum region. Such a transition has to 
be performed with the due account of the Tauberian theorem. This comment relates 
to the recent ALPHA collaboration results on the asymptotic behavior of the QCD 
effective coupling obtained by lattice simulation. 
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1 Invariant coupling in different representations 



Current practice of the QFT calculations, as a rule, employs expression of observables and 
of some other intermediate renormalizat ion-group (RG) invariant, or covariant, objects in 
terms of the invariant coupling function (also referred to as the "invariant, or running coupling 
constant" 1 ). This invariant coupling a(Q 2 ) is a real function of a real positive argument 
Q 2 = Q 2 — Ql , momentum transfer squared. 

The notion of invariant coupling (IC) has initially been introduced 2 — see, e.g., refs.^Ej 
on renormalization group — in terms of a product of real constants Zi entering into finite 
Dyson renormalization transformations 

Di(...,a) -> £>■(..., a') = z^ 1 Di(...,z 3 a>) ; a^a' = z 3 a, (1) 

with particle propagators £)*(..., a) = Di(Q,m, /x; a) taken in the energy-momentum rep- 
resentation. Here, the IC is expressed in terms of scalar (i.e.,, Lorentz-invariant) QFT 
amplitudes taken also in the momentum representation. For instance, in QED one has 

a{Q )=ad ^—,—- a j 

with d , the transverse photon propagator amplitude which is real in the Euclidean domain. 
In QCD, IC a s (Q 2 ) is usually defined in a similar way as a product of the expansion 
parameter a s with an appropriate scalar vertex and propagator amplitudes, see, e.g., [3| 
and references therein. 

This IC is commonly used for parameterization of the renormalization-invariant quan- 
tities, in particular, of observables. To this goal, one should take an invariant quantity of 
interest in the appropriate representation. In particular, in terms of ct s (Q 2 ) there could be 
presented only an object taken in the momentum representation. Moreover, this object to 
be expressible in terms of the real function a s (Q 2 ) should be real itself. In the QFT case, 
this corresponds to the space-like Euclidean domain with Q 2 > . 

Meanwhile, some observables, like effective cross-sections, are functions of another Lorentz- 
invariant argument s , the center of mass energy squared. As it is well known, the polarization 
operator H(Q 2 ) , being represented in the form of the Cauchy type integral, provides us with 
a technical means for the integral connection of invariant functions of the Q 2 and s ar- 
guments. Its logarithmic derivative, known as the Adler function, is defined via integral 

close to the Kallen-Lehmann spectral representation for H(Q 2 ) . Here, R(s) being the imag- 
inary part of the polarization operator, due to the optical theorem, is proportional to the 
total cross-section which is an observable quantity. 

x In view of semantic absurdity of the last term, we use the expression "invariant coupling function" or 
"invariant coupling" . 

2 Under the name of "invariant charge" (of electron), the natural one in the QED context. 
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Further on, we shall take relation (J2J) as the base for definition of the integral operation 
D transforming a real function M(s) of the positive (time-like) argument into another real 
function E{Q 2 ) of the positive real (space-like) argument Q 2 — see below eq.(j7J). 

More generally, we consider a linear integral transformation with kernel K(x) , depending 
on one argument 



that comprises various kinds of one-dimensional Fourier transformations (usual, Laplace, 
cosine-Fourier, sine-Fourier and some others, e.g., with a kernel in the form of the Bessel 
function that could be derived from the four-dimensional Fourier one) and, with K(z) = 
(1 + z 2 ) 2 , just mentioned the "Adler transformation" D defined below by eq.(|7|). 

In Section 2, we start with the discussion of the interrelation between the RG invariance 
and integral transformation. Our first explicit example taken from the so-called Analytic 
Perturbation Theory (APT) concerns a transition from the Minkowskian energy-momentum 
domain to the corresponding Euclidean one. We continue with the transition from the "dis- 
tance" (i.e., , three-dimensional space) representation to the common momentum (transfer) 
one using the one-dimensional sine-Fourier transformation. 

Further on, in Section 3, we shall pay special attention to the relation between a long- 
range distance and infrared momentum asymptotic behaviors in QCD, and comment on the 
interpretation of some recent results of the ALPHA collaboration. 

2 Integral transformations 

2.1 RG invariance through integral transformation 

First of all, note that transformation L, eq.(3), is compatible with the RG transformation. 
For example, in the massless case, let some initial function f(x/fi;a) be invariant with 
respect to the RG transformation 



Then, its integral image F(yfi; a) = h[f](y) will be invariant with respect to the same RG 
transformation (@J) — and vice versa — with the same transformation function a tr which 
satisfies the functional equation 



that follows from the group composition law RqR t = Rq t . Note here, that eq. (jSJ) implies the 
"canonical normalization condition" for the transformation function a tr 




(3) 




(4) 




(5) 




(6) 



that corresponds to the identity transformation R\ . 
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Relation (J2J) is commonly used to define the Adler function D . However, we shall treat 
it instead as a definition of the integral operation D transforming a function M(s) of the 
positive real argument s into another function E(Q 2 ) of the positive real argument Q 2 



r 

M(s)^B[M}(Q 2 ) = Q 2 / 



oo 



ds 



(s + Q 2 ) 2 



M(s) = E{Q 2 ) 



(7) 



with the reverse operation R = [D] _1 that can be expressed in terms of the contour integral; 
for more details about these operations, see, e.g., ref.^Q Ej. 

The transformations D and R preserve the RG invariance; for instance, they transform 
a renormalization invariant function M(s/n 2 ;a) into another RG invariant E(Q 2 / fi 2 ; a) 
and vice versa. The first of these invariants M should be representable as a function of 
an adequate "Minkowskian invariant coupling" a(s) , while the other E — via the common 
invariant coupling a(Q 2 ) . In particular, this means that starting with the usual QCD cou- 
pling 6t s (Q 2 ) , by the operation R it is possible to define 3 the QCD effective coupling in 
the time-like Minkowskian region — R[«s]( s ) m terms of a s (Q 2 ) , a QCD invariant 

coupling initially introduced in the Euclidean domain. 

Another example is related to the physical amplitude A depending on a couple of Lorentz 
invariant kinematic arguments s and t . Let it be renormalization invariant 



Here, there are two possibilities. The first one deals with an integral transformation with 
respect to one of the two kinematic arguments like an eikonal transformation F(s, t) — > 
$(s, b) . A new amplitude $ has the RG transformation property similar to that of F . Like 
in the general case of several kinematic arguments, the RG transformation properties of F 
and $ are not very useful. 

However, this situation changes for an integral transformation of the second type involv- 
ing some function of ratio of both the variables (p(s/t) . Such a transformation 



projects the initial function of two (kinematic) arguments onto a (set of) function(s) A k (s) 
of one argument which is(are) RG invariant. The transition to partial waves and to moments 
of structure functions provides us with examples. 

One more example of an integral transformation of the "one-argument" function compat- 
ible with the RG invariance is given by the transformation of the Fourier type which follows 
the general linear form, eq.flSJ). It relates the function f(r) of the space-time Lorentz in- 
variant argument r = \/v 2 — t 2 with the function F(Q) of an energy-momentum invariant 
argument Q = Q 2 = Q 2 - Ql . 

3 However, in practice, the nonphysical singularity in perturbative a s (Q 2 ) prevents one from a straight- 
forward performance of the integration procedure. In more detail, this issue has been discussed in [5]. 



A(s, t; ft, a) = A(s, t; fjf, a') . 




■oo 



A(s, t) K (k) (s/t) d<p(s/t) ; k = 0, 1, 2, ... . 
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As it is well known, the Dyson transformation (fl|). being "a self-similarity property of the 
Schwinger-Dyson equations" , can be, on equal footing, considered in both the (energy- 
momentum — like in eq.Q — and (time-) space representation, that is for functions like 

D(...,a) = D(x,m, n;ot) ; x = {r, £}. 

In the QED context, this latter picture 4 was used by Dirac while discussing polarization 
of vacuum in terms of effective charge of electron first introduced by himjZj as a space 
distribution 

/ \ f a o , r 1 e n 1 h , . 

e D {r) = e U - — In — \ ; a = — ~ 1Q7nQ „ n r e = 8 

[ 3n r e ) 4n 137.0360 m e c 

of the electric charge around the point "bare" electron. 

More recently, an analogous object, the QCD effective coupling g 2 (L) , which is a function 
of the spatial size L of a lattice, has been introduced|H] and used[§] [12] by ALPHA 
collaboration. 

Some other examples of the relation between integral transformations and renormaliza- 
tion group symmetries in problems of classical and mathematical physics can be found in a 
fresh review paper |T3| (see, especially, Example 2 on page 358 and references therein). 

The issue of correlation between the RG formulations in different pictures has two aspects. 
The first one deals with the basic notions and objects of RG transformations. The second 
one is that of these objects correlation via an appropriate integral transformation. 

According to the terminology formulated in Refs.|14[ll5j. an invariant function satisfying 
the functional equation (j^J) and canonical normalization condition should be named effective 
coupling (EC). To the case of an invariant function with a more general normalization 

a) = N(a) ^ a 

there corresponds a term invariant coupling (IC). Contrary to EC, it can not be used as 
a function transforming a coupling constant in (j3J). Generally, integral transformation (jSJ) 
maps an EC a(x) onto some other RG-invariant function A(y) which is IC rather than EC. 

A few questions are in order: 

a) How to define a and its integral image A explicitly for the given QFT model? 

b) How to relate them? 

c) Which of them can be used as the transformation function a tr ? 
"Immediate" answers — 

a] Use a common algorithm with perturbative beta-functions to define a and A , 

b] Relate them by an appropriate transformation (JHJ) 
— turn out to be incompatible with the each other. 

4 We shall refer to it as to the "distance representation" . 
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2.2 From the cm. energy to the momentum picture 

This situation can be illuminated by the known answers to the same questions for the 
integral transformation Q. Let us write it down for the particular case of the coupling 
function transformation 



a(s) 



N(g 2 ) = Q 2 / 

Jo 



ds 



[s + Q 



2\2 



at(s) 



a(Q 2 ). 



(9) 



It turns out to be impossible to use here a common perturbative QCD coupling a s {Q 2 ) 
as ot{Q 2 ) because of its unphysical singularities, like the Landau pole at Q 2 = A 2 , which 
contradicts the integral expression (J0J). The latter implies that function a(Q 2 ) should be free 
of any singularities outside a cut > Q 2 > — oo . The same is true for the candidature of 
a s (Q 2 ) for the role of function a(s) . One of the possible solutions that was proposed in the 
so-called Analytic Perturbation Theory (APT) (see, refs. [T7J consists in using of a s {Q 2 ) 
only as a prototype for both the u(Q 2 ) and a(s) which in the one-loop case with 



a 



(i) 



(Q 2 



/3 hi(Q 2 /A 2 ) 



A 



\x exp 



turn out to be 



als) 



arccos 



7T- 



and 



«(Q 2 ) = W 
Po 



1 
L 



L s >0 



A 2 



arctan(7r/L s 
n(3 



In 



A 2 



Q 2 - A 2 



T 1 ^ 2 

A 2 ' 



(10) 



11' 



(12) 



Both these two functions, corresponding to (fTUj) in the weak coupling limit, are ghost-free 
monotonous functions related by the transformation (JHJ) <y{Q 2 ) = ^[o]{Q 2 ) and its reverse. 

As it has been noticed in Ref.jH], transitions from expression (fTQj) for a { s\Q 2 ) to a(s) 
and Oi{Q 2 ) can be represented as a consequence of a transition from the usual QCD coupling 



constant a s to the new ones 



O'.s- 



and 



arccos 



a E 



O'.s- 



arctan(7r/3o a s 



a.+i-(l-e 1 /A«.)- 1 
Po 



Hence, a transition from the Minkowskian coupling a(s) to the Euclidean one ot{Q 2 ) is 
equivalent to that one induced by the following coupling constant transformation : 



oc M — ► "£(«*/) 



— tan(7r/5 a A /) + 3- • 7- 



3 7rcot(7r/3 OM) 



(13) 



at < «m < 1/Po and 0<a£<l//3 - 

In turn, this implies that the integral transformation generally changes the normal- 
ization of coupling functions and, in particular, maps the EC onto IC. 
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2.3 Momentum and distance representations 

Turn now to the Fourier transformation. It will be convenient to use the modified sine- 
Fourier transformation 5 of the RG-compatible form 

2 r°° dr 

F sin/r [f](Q) = - / — sin(Qr) /(r) = F(Q) (14) 
71 Jo r 



and its reverse 



POO 

FTMF](r) = r / dQ sin(Qr) F(Q) = /(r) . 
./o 



Now, rewrite expression (jHJ) squared in a more contemporary and general notation 



2 \ z' 2 

r \ \ a„ , r 



ag w ,a v =a n 1 "^ ln ^i (15) 

(r M being a reference distance) with the correspondence relation ct^ = ao at r M = r e . This 
distance-representation perturbative EC a p ^(p 2 ,a) can be connected by the transformation 

a(Q 2 /p 2 ,a) = F sin/r [ag](Q 2 ) 

with a more common perturbative QED coupling in the momentum representation 

a P i — ,a M =a M <l + — In— f with p = — (16) 



and eg = e c = 0.5614 . Expression (JTfij) can also be obtained from (fTB^) by the argument 
substitution 

r -> c E /Q , (17) 

that is close to the quantum-mechanical correspondence relation r — > 1/Q being slightly 
modified by changing the Q scale. By the way, the positive effect of the cf factor is clearly 
seen in our Figure 2 — see below. 

Quite analogously with the presented example, both the functions a RG and a RG , being 
explicitly defined (in the familiar form of the geometric progression sum) by the common 
RG perturbation-based procedure with (fT5]l and (fT^|) as inputs, are not related by eq. (fT4^) . 
Instead, they are connected by a more involved relation a RG (Q) = ^ {F sin [a^ G ] (Q)} ■ 

Moreover, generally, Fourier transformation maps an effective coupling onto some in- 
variant coupling that cannot be used as a tr in (j3J). For example, if one starts with the 
perturbative distance-representation EC of the form 



2 



a% (p 2 ,a) = a - (3 Q (a) hip 2 + a {J^-} ln 2 p 2 ; p 2 = j^; a = 



3 that is a particular case of J3J with K(z) — 2sin(z)/7rz and follows from the usual 3-dimensional one 
$(Q) - (2tt)- 2 / driP{r)e lQr with F(Q) = Q 2 Tp(Q) , f(r) = r ^(r) . 
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then, as a result of the sine- Fourier transformation (|14|1. one arrives at 



/ Of \ 2 

a pt (g 2 ,a) = a + /3(a) lng 2 + a (^— J {ln 2 g 2 + A 2 } (19) 

with 



Q 2 _ ( Qth\ a _vr 2 



Hence, a(l,a) = a + a(a/37r) 2 A 2 = a' 7^ a and a(g 2 ,a) is an invariant coupling 6 . By 
transition to another coupling constant a — > a' it can be transformed (see Section 3 in 
into the EC a'(q 2 ,a') . The transformation is relevant to the ty(a) dependence. This could 
be important at the strong coupling case in the IR region. 

Note also that the logarithmic terms in both ()18|) and (|T9"j) yield the sums of geometric 
progression with arguments of logarithms related by eq.([17j). 



3 Long— range and infrared asymptotics 
3.1 Tauberian theorem 

Turn now to the particular issue of correlation between asymptotic behaviors of the functions 
/ and / related by the Fourier transformation. 

This correlation is popular in quantum mechanics where, quite often, one uses the so- 
called "quantum-mechanical correspondence relation" 

r^l/Q, (QMC) 

which in the IR case is equivalent to 

F{Q)^f{Q- 1 ) as Q^O. (20) 

Heuristically, this last feature could be simply understood by a change of the integration 
variable r — > x = rQ in the general linear transformation (JHJ) 

*■(«) = /"£*(*>/(!). (3a) 

However, for a more rigorous derivation of ()20|) one needs to specify some asymptotic 
property of the function f(r) as r — > 00 . In short, this can be formulated as the Tauberian 
theorem 7 : ( Here, the symbol " ~ " means "behaves like".) 

6 Meanwhile, this feature is not essential (like for transformation (|13f) ) in the weak coupling case at the 
one- and two-loop levels. 

7 Originally, under the name of Tauberian theorems one implied statements concerning the relation be- 
tween summability and convergence of series. More recently, in the middle of XX century, this term started to 
be used in the context of asymptotic properties of integral transformations. Here, we give only a crude out- 
line of this important theorem for the Fourier transformation, the sketch that is sufficient for our application. 
For a more complete and rigorous exposition of this matter the reader is referred to refs. [261 12 7j . 
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If function f(r) asymptotically satisfies "the separability condition" 

f(kr) ~ Ccj)(r) g(k) as k -> oo with C ^ , (5) 
then, under some additional conditions, its Fourier image obeys the property 

F(Q)~Q{1/Q) as Q^O, (T) 

— that, with some reservation, follows from eq.(3a). 

Now, for a definite class of functions g(k) entering into the condition (S), e.g., of power 
and/or logarithmic type 

f( r ) ~ f?( r ) ~ r ^ (lnr) 7 (In lnr) 5 . . . as r — > oo , (21) 

it is possible to obtain from (T) the correspondence rule (|2()j) . 

3.2 Illustrations 

Consider few examples given in the Table. 

Table 





f(r) 


F(Q)=¥ sin/r [f] 




(20) m 


{*} 


rF(r) 


2/(Q)/ttQ 


r — > oo 


g^o 


{1} 


r v \ o < v < l 


A sin (-)r(l + Ber- 


+ 


+ 


{2} 


lnr 


ing- 1 -C; C = 0.5772 


+ 


+ 


{3} 


(lnr) n ; n > 2 


(lng-i-Cr + A, 


+ 


+ 


{4} 


r 2 
r 2 +o? 


e" aQ ; a > o 


+ 


+ 


{5} 


e~ ar ; a > o 


- arctan f — ) 

it \ a 1 







In two right columns we mark the correspondence of the function /(r) at r — > oo to the 
Tauberian condition (S), as well as the fulfillment of condition (QMC) only in the IR region 
in the form (|2()|). Note, that symbols "+" and "-" turn out to be completely correlated. 

As it follows from the Table, in accordance with eq. (j21j) . the class of admissible functions 
is rather narrow. For instance, it does not contain trigonometric functions and exponentials 
~ e ar . Only the first-line expression F^ij corresponds (up to a factor!) to the (QMC) rule. 
The second line's one _F{ 2 } = ln(g _1 ) — C, (with C being the Euler constant) satisfying 
the Tauberian condition as Q — > , generally differs from /{2}(l/g) , especially in the region 
Q ~ 1 . The same is true for the line {3}; here the constant A 2 had been introduced before 
in eq. ffTHjl. 

The next line {4} provides us with an example that satisfies the Tauberian condition (S), 
but severely violates the (QMC) rule for the whole function F^y. It is instructive to compare 
the behavior of the Fourier image F^{Q) = exp(-aQ) with the "QMC-substituted" initial 
function /{4}(l/g) = (1 + a 2 Q 2 )~ 1 . As it can be seen in Figure 1, the relative error being at 
aQ > 0.2 of an order of 15%, quickly increases and reaches the level of 30% at aQ ~ 1 . 
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-2-10 1 2 



Figure 1: Behavior of the "QMC-substituted" function /{4}(1/Q) = (1 + a 2 Q 2 ) 1 and of its sine-Fourier 
transform F^y(Q) = exp(— aQ) . 



Only last example /{ 5 } does not satisfy the Tauberian condition, while the corresponding 
limiting values /{ 5 }(oo) = F{ 5 }(0) = coincide. 

To give some example more close to the realistic QCD case, consider the sine-Fourier 
transformation for the class of functions f(Q) that satisfy the Kallen-Lehmann spectral 
representation 

2 [°°dQ P dcrpja) 

/(?) F ( r ) = - / "77" sra(Qr) / (22) 
k Jo Q Jo v + Q 

with the density p(er) . Changing the order of integration and performing the integration over 
Q with the help of line {4} from the Table we arrive at 

F(r) = J™ P -^- (l - e-V*) . (23) 

In the case of positive spectral density, this expression defines a monotonically rising 
function of r with possible singularity as r — > oo . 

As an explicit example, we consider the analyticized invariant QCD coupling a(Q 2 ) , 
eq.flHJ with the spectral density Papt^) ~ [(lncr) 2 +7r 2 ] _1 taken from the one-loop APT [IS]. 
For oapt = A)«apt this yields 



f°° da(l - e" r v^ T . . r . . f°° e~ r ^da 

aAPT (r) = / -r^ ^ 57 = 1 - I r ; 7(r = / r= =r . 24 

7 aKln^s + vr 2 ] 7 a[(ln a) 2 + vr 2 ] 

The r.h.s. integral J(r) resembles the Ramanujan ones|19j 

/■oo e ~ TX dr r°° p~ t ^° do 



x[(lnx) 2 + 7r 2 ] Jo a[(lna) 2 + 47r 

expressible in terms of the special transcendental function u(x) . We shall use this proximity 
for analysis of the J(r) behavior as r — > . The function v asymptotics is well known [TH] . 
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One has 



The difference 



R(r) 



lnr 2 



In 2 r 2 ' 



A(r) = I(r) - 



R(r) 



, da 
3tt z / — 



as r — > . 



(25) 



2 7 a [(l n( T) 2 + 7r 2 ][(lncT) 2 + 47r 2 ] 

being positive, vanishes A(oo) = at infinity. At the same time, A(0) = 1/4. Hence 

1 



a A pT(r 



In r 2 



as r — ► and aAPT(oo) = 1 . 



(26) 



It is instructive to compare «apt with the monotonous initial function 



a an (Q 2 ) = P a(Q 2 

Oan(0) 

in the coordinates r = 1/Q . 

1.0 



1 1 

+ 



o 



w 2 + Q 2 ln(Q 2 ) 1-Q 2 
1 



InQ 2 



as Q — > oo 



(27) 



a D d /q) 




iQioQ/A 



Figure 2: Behavior of the APT-analyticized a an (Q) and of its Fourier transform od(1/Q) with the 
argument modified by the quantum-mechanical rule" (QMC) and of its Fourier image od(cf/Q) (dotted 
line) with the scale modified by (17). 

Such a comparison (see Fig. 2) reveals a puzzling similarity of both the functions od(1 / Q) 
and a an (Q 2 ) , that could be essentially improved by changing a^l/Q) for a D (cp/Q)- In 
particular, the relative error in the region 4 < Q/A < 70 is reduced from 30% to 10%. 
This produces an impression that one has a strong argument for supporting the use of 
the correspondence rules (QMC) and (fTTj) for the qualitative estimate of Fourier image 



11 



of functions with an "adequate" - - to eq. (}2*T]) - - asymptotic behavior. Unhappily, this 
impression is erroneous. For example, the analysis of the family of functions 



with adequate asymptotics reveals 8 that similarity like in Figure 2 represents a "puzzling 
exclusion" of usual case, which, typically, is far even from very approximate similarity like 
the one in Figure 1. 

Hence, the (QMC) rule, generally, does not provide us with some reasonable means for 
qualitative estimation of the Fourier image behavior in the region far from the asymptotic 
one even for functions satisfying the Tauberian condition. 



3.3 The "Schrodinger functional QCD coupling" 

To analyze by nonperturbative means the infrared behavior of QCD, the "ALPHA collabo- 
ration" uses the functional integral approach (both in the quenched QCD version and with 
two massless flavors). It works with the Schrodinger functional (SF) defined in the Eu- 
clidean space-time manifold in a specific way: all three space dimensions are subject to 
periodic boundary conditions, while the "time" one is singled out — the gauge field values 
on the "upper" and "bottom" lids differ by a phase factor with the parameter 7] . Then, the 
renormalized coupling osf is defined via the derivative T' = dT/dn of the effective action 
T = a~ l Y Q + T 1 + a s T 2 + ■ ■ ■ as (cf. eq.(8.3) in ref. [TU]) a function in the distance represen- 
tation a SF (L) = T' /T' with L being the spatial size of the above-mentioned manifold. 

Quite fresh results reveal the steep rise of the SF asp(L) coupling with L in the region 
«sf — 1 • Here, the analytic fit ^2] to the numerically calculated behavior of cvsf has an 
exponential form 

a SF (L)^e mL (28) 

with m ~ 2.3/L max and L max (a reference distance in the region of sufficiently weak 
coupling) indirectly defined via the condition asF(-^max) = 0.275. For a discussion of the 
momentum-transfer QCD behavior, in the papers of ALPHA collaboration the "quantum- 
mechanical correspondence" rule is used in the form 

L -> l/n {QMC 1 ) 

with fi practically treated as a reference momentum transfer value. This way of transition 
from the distance to the momentum transfer picture is equivalent to the one discussed above 
in Section 3.1. It works quite well in the UV region as far as ~ In -1 Q 2 behavior is compatible 
with the Tauberian theorem. Nevertheless, in some precise consideration, like in numerical 
relation between various scales [9 j 9 , one should take into account the modification of scale 
according to eq. lfTTJ) . 

8 Private communication by Dr. A.V. Ncsterenko. 

9 This important paper was brought to our attention by Dr. U.Wolff. His friendly assistance is gratefully 
acknowledged. 
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At the same time, the steep rising function of the type (J28j) does not fit to the (S) 
condition of this theorem. Due to this the IR (i.e., as Q 2 —>■ ) behavior needs a more 
elaborate analysis. 

3.4 On the representation of the ALPHA results in the IR region 

Generally, there are two different standpoints for discussing this transition for numerical 
simulation results obtained on a finite lattice of size L . 

The first, simple-minded one treats the obtained numerical results and their analytic 
description (|28|) as an approximation to some limit with the continuous function a SF (L) 
that is subject to an integral Fourier transformation. On the other hand, as far as a lattice 
boundary condition is a periodical one, it is possible to use a Fourier series instead of the 
Fourier integral. 

In what follows, we shall try to discuss the possible Fourier image of asp(L) having in 
mind both the possibilities. 

To start with the integral Fourier transformation, we represent the ALPHA "distance 
running coupling" q;sf by the sum of two terms: 



with «pt , a perturbation contribution, and «al , an essentially nonperturbative part. 

As far as results of all nonperturbative calculations — by both the numerical lattice 
simulation ones and solving the Schwinger-Dyson equations (SDE) — reveal no traces of 
unphysical singularities, we, on the one hand, change the first perturbative term apT for 
some regular expression optf (like the "freezed" one (|24JI in the Analytic Perturbation 
Theory or some other smooth ones like those emerging from the "effective massive glueball 
model" [201 or from the SDE solving EU). 

On the other hand, we approximate the second, essentially nonperturbative, term by 
expression of type 



It is close to (J28j) and admits further explicit integration. Here, k > 1 is small integer, and 
the mL m ~ 2.0 and Ck paremeter values follow from results of numerical lattice simulations. 
For example, G\ ~ 0.013 . 



For the Fourier image of regularized ttpT— r we assume, qualitatively, the momentum 
transfer behavior related by the (QMC) rule 



«sf(£) = «pt(£) + cxkl{L) 




(29) 



That is 



<*sf( l ) = «ptf(£) + a k {L) . 



(30) 



ort-r(Q) = F sin [apTF](<3) = «pt-r(1/Q) 



(31) 



with a finite IR limit like in (|27|) 



o;pt-r(0) = C ; < C < oo 



(32) 
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which is supported by the second illustration of Section I3~2l see Figure 2. 

To perform transformation (fT4*j) for aifc(L) , one needs some regularization. To this goal, 
we shall insert the cutoff factor 10 

exp{(m + n)6(L-CL m )^L m -L)}; £ = 2.5-^3, p > 

into the integrand of the r.h.s in eq. (JHJ). Here, parameters £ and fi , regulate "the place of 
switching-on" and the intensity of the cutoff. 

Then, the "ALPHA coupling function" in the momentum picture can be represented as 
a sum of three terms 

«Sf(Q) = "ptf(Q) + ak,rcg(Q) + "k.singW; n) (33) 

with the first one being defined by (}3"Tj) and two others can be obtained from expressions 
1(25)1 taken at k = 1 : 

ai,reg(Q) = y 1 / ^sin(QL)e mL = -1 [/i(Q,£L m ) - h(Q,L m )} ; (34) 

«WQ;/i) = e m ^-^ / dLsin(QL)e^ Lm - L ) = -J-^(Q, £L m ; M ) (35) 
by appropriate differentiation with respect to m . Here, 

/i(Q L) = e mL msin (^ L ) -Qcos(QL) 



m 2 + Q 2 
mL Qcos(QL) + Aisin(QL) 



In the case of the periodical function with a period related to £L m the transformation 
results in the Fourier series rather than the Fourier integral and the third term in (}3*3"|) is 
absent. 

The most important qualitative feature of the functions fk(Q,L) and <fk(Q, L; //) is 
their IR behavior. All the functions generally tend to zero linearly with Q — > . The only 
exception is the case of /x = when one has a power singularity (pk(Q — > 0, L; 0) — > 1/Q . 

By combining this result with (|32|) we conclude that «sf(Q) * n the ^ region can have 
a finite limit or the first order pole (the latter - - only for the case of integral Fourier 
transformation), contrary to the exponential growth ~ exp (m/Q) that could be anticipated 
from some results of the ALPHA collaboration — we mean, e.g., transition from Fir. 4 to 
Fig. 3 in Ref.|E|. 

On a more general ground, one can argue that the exponential growth of the QCD 
coupling with L is the utmost steep possible one. In particular, for the In ct(L) ~ L v , v > 1 
regime there is no known mathematical means for defining a Fourier transformation. Vice 
versa, for the ~ exp{m/Q} IR asymptotic behavior it is impossible to construct any Fourier 
transformation and "return" to the L-picture. 

10 Strictly speaking, this exponential regularization corresponds to the generalized Fourier transformation 
— see, e.g., Section 1.3. in ref.|22). 
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4 Discussion 



Our first observation concerns the issue of unphysical singularities in the QCD effective 
coupling and related QFT quantities like the propagator amplitudes. These singularities 
are inconsistent with integral transformations. In particular, the well-known unphysical 
singularity in the region Q ~ A ~ 300 — 400 GeV prevents the common QCD perturbative 
effective coupling a s (Q 2 ) from a straightforward performing of the integration procedure 
necessary for transition to the distance representation by the Fourier transformation. This 
gives us an additional theoretical argument against of existence of unphysical singularities 
in the physically reasonable theories like QED and QCD. Remind here, that all the QCD 
nonperturbative calculations — by both the numerical lattice simulation and the SDE solving 
— reveal no traces of unphysical singularities. 

The second result based on the analysis of the Fourier transformation deals with the 
quantum-mechanical correspondence rule 

r -> 1/Q (QMC) 

relating the asymptotic behavior of a function f(r) as r — > oo and of its Fourier transform 
F(Q) as Q — > . It has been confirmed that this rule, being a reasonable guide for some 
class of asymptotics (the power and logarithmic type), has its rigid limits of applicability. 

First, even for the function with admissible asymptotic behavior, in the region not very 
close to the singularity the (QMC) rule yields only a qualitative correspondence, as it follows 
from our Figure 1. 

Second, it is not valid at all for wide class of asymptotic behaviors violating the so-called 
Tauberian conditions |27j. like the exponential ones. 

In particular, the exponentially rising long-range behavior of QCD coupling in the dis- 
tance representation 

a SF (L) ~ e mL 

observed by ALPHA collaboration on the basis of the lattice simulation of the Schrodinger 
functional, according to our analysis can correspond in the momentum (transfer) picture to 
the 

a) finite or b) "slightly singular" ~ 1 /Q (ALPHA-IR) 

IR asymptotics. 

This means that these long-distance results, being properly translated to the IR region, 
qualitatively, will not be so far from the results of other groups 11 that perform lattice simula- 
tion calculations (partially supported by solution of appropriate truncated Schwinger-Dyson 
equations) for the functional integral defined in the momentum representation. 

From the physical point of view, in our opinion, there are at least two issues that should 
be mentioned in connection with the QCD infrared asymptotic behavior. 

First, we have to remember that the region of Q < 500 MeV, physically, corresponds 
to distances r > 10 _13 cm , that is to the hadronic scales. Here, all the QCD notions, like 

11 See, e.g.. refs. |21l 1231 and short discussion of their difference in |15| . 
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gluonic and quark propagators, seem to be meaningless. Even in a more moderate region 
0.5 MeV < Q < 3 GeV of strong QCD interaction there arises a question of physical meaning 
of nonperturbative QCD functions, including the effective coupling one. 

The mentioned above results for the effective QCD coupling obtained [2U EHl 121] by a 
numerical lattice simulation of the path QCD integral in the momentum representation are 
formulated by the "common QFT language" using the vertices with fixed dynamics, like, 
e.g., in [231213 for the QCD model with two massive quarks. There, an invariant coupling 
g{Q 2 ) is defined on the basis of the gluon-quark vertex T q _ g i(0; Q 2 , Q 2 ) in the particular 
MOM scheme with gluon momentum equal to zero. The invariant coupling function thus 
defined suffers from the usual drawback of MOM schemes — the gauge dependence. Never- 
theless, it can be, in principle, considered as directly corresponding to some definite physical 
situation (being incorporated into a series for some observable with gauge-dependent coef- 
ficients) . 

Second, it seems to be reasonable to relate the IR behavior of the QCD functions with 
the confinement phenomenon. Here, it is possible to appeal to the so-called Kugo-Ojima 
condition 28] that, physically, corresponds to the absence of "open colour" in the asymptotic 
states. In the QCD language, this yields the vanishing of the gluon and quarks fields renor- 
malization constants, that, in turn, is equivalent to the zero IR limit of corresponding QCD 
propagators. Such a behavior has been observed in the most of the "lattice-simulation QCD 
papers" j2H 123 121] ■ Quite recently it has been supported by Orsay group [29] on the basis of 
an instanton liquid model and, in a sense, by analysisj^H] of the r decay data. 

In this context, the (ALPHA-IR,a) IR behavior seems to be a quite reasonable possibility 
to correlate all above-mentioned lattice simulation results with the each others and with the 
physics of confinement. 
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